Low Momentum Nucleon-Nucleon Interactions and Shell-Model Calculations 
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In the last few years, the low-momentum nucleon-nucleon (NN) interaction Vi ow _k derived from 
free-space NN potentials has been successfully used in shell-model calculations. Viow-k is a smooth 
potential which preserves the deuteron binding energy as well as the half-on-shell T-matrix of the 
original NN potential up to a momentum cutoff A. In this paper we put to the test a new low- 
momentum AW potential derived from chiral perturbation theory at next-to-next-to-next-to-leading 
order with a sharp low-momentum cutoff at 2.1 fm _ . Shell-model calculations for the oxygen 
isotopes using effective hamiltonians derived from both types of low-momentum potential are per- 
formed. We find that the two potentials show the same perturbative behavior and yield very similar 
results. 



PACS numbers: 21.30.Fe, 21.60.Cs,27.20.+n,27.30.+t 



I. INTRODUCTION 

A challenging goal of nuclear structure is to perform 
shell-model calculations with single-particle (SP) ener- 
gies and residual two-body interaction both derived from 
a realistic nucleon-nucleon (NN) potential Vnn- To 
tackle this problem, a well-established framework is the 
time-dependent degenerate linked-diagram perturbation 
theory as formulated by Kuo, Lee and Ratcliff [l], Q], 
which enables the derivation of a shell-model effective 
hamiltonian H c q starting from Vnn- 

As is well known, a main feature of Vnn 1S the pres- 
ence of a built-in strong short-range repulsion dealing 
with high-momentum components of the potential. This 
hinders an order-by-order perturbative calculation of H c g 
in terms of Vnn, as the matrix elements of the latter are 
generally very large. However, in nuclear physics there 
is a natural separation of energy scales, that can be used 
to formulate an advantageous theoretical approach for 
Vnn- As a matter of fact, the characteristic quantum 
chromodynamics (QCD) energy scale is Mqcd ~ 1 GeV, 
while for nuclear systems we have M nuc ~ 100 MeV Q. 

This consideration was at the origin of the seminal 
work of Weinberg 0, who introduced into nuclear 
physics the method of effective field theory (EFT) to 
study the S'-matrix for a process involving arbitrary num- 
bers of low-momentum pions and nucleons. This ap- 
proach is based on the known symmetries of QCD and 
parametrizes the unknown dynamical details introducing 
a number of constants to be determined. 

Since then much work has been carried out on this 
subject (see for instance 0, 0, 0, H S leading to 

the construction of Vn n 's based on chiral perturbation 
theor y t hat are able to reproduce accurately the NN 
data [mill]. However, also these potentials cannot be 
used in a perturbative nuclear structure calculation. 



Inspired by EFT a new approach based on the renor- 
malization group (RG) has been recently introduced to 
derive a low-momentum AW interaction V[ ow _k [HI, [HI, 
Il5j |. The starting point in the construction of Vj nw -k 
is a realistic model for Vnn such as the CD-Bonn (lq . 
Nijmegen [l7j], Argonne V18 or N 3 LO [U poten- 
tials. A cutoff momentum A that separates fast and slow 
modes is then introduced and from the original Vnn an 
effective potential, satisfying a decoupling condition be- 
tween the low-and high-momentum spaces, is derived by 
integrating out the high-momentum components. The 
main result is that Vj ow -k is a smooth potential which 
preserves exactly the onshell properties of the original 
Vn n i and is suitable to be used directly in nuclear struc- 
ture calculations. In the past few years, Vi ow -k has been 
fruitfully employed in microscopic calculations within dif- 
ferent per turbative frameworks such as the realistic shell 
model [SIM El iHlMlllithe Goldstone expansion for 
doubly closed-shell nuclei [25L l2fl HtJ , and the Hartree- 
Fock theory for nuclear matter calculations [M [M [3(| • 

The success of Vj ow -k suggests that there may be a 
way to construct a low-momentum NN potential which 
is more deeply rooted in EFT and can be used in a 
perturbative approach. Namely, instead of taking the 
detour through a NN potential with high-momentum 
components, one may as well construct a low-momentum 
potential from scratch using chiral perturbation theory. 
We have constructed such a potential at next-to-next-to- 
next-to-leading order (N 3 LO) using a sharp cutoff at 2.1 
fin -1 . This potential reproduces the NN phase shifts up 
to 200 MeV laboratory energy and the deuteron binding 
energy. While Vj w-k allows only for a numerical rep- 
resentation, the new low-momentum potential (dubbed 
N 3 LOW) is given in analytic form. Moreover, the low- 
energy constants are explicitly known so that the chi- 
ral three-nucleon forces consistent with N 3 LOW can be 



2 



properly denned. 

In order to investigate the perturbative properties of 
this new chiral low-momentum potential, we perform 
shell-model calculations for even oxygen isotopes. We 
employ two different effective hamiltonians: one is based 
on the Mow-k derived numerically from a 'hard' N 3 LO 
potential [ll| and the other on N 3 LOW. 

The paper is organized as follows. In Sec. II we give a 
short description of how Vi ow -k is derived as well as an 
outline of the construction of the N 3 LOW potential. In 
Sec. Ill a summary of the derivation of the shell-model 
H e f[ is presented with some details of our calculations. 
In Sec. IV we present and discuss our results. Some 
concluding remarks are given in Sec. V. 



II. LOW-MOMENTUM NUCLEON-NUCLEON 
POTENTIALS 



Q is an irreducible vertex function whose intermediate 
states are all beyond A and Q' is obtained by removing 
from Q its terms first order in the interaction Vnn- In 
addition to the preservation of the half-on-shell T matrix, 
which implies preservation of the phase shifts, this VJ ow -k 
preserves the deuteron binding energy, since eigenvalues 
are preserved by the KLR effective interaction. For any 
value of A, the low-momentum potential of Eq. ^ can 
be calculated very accurately using iteration methods. 
Our calculation of V[ ow _k is performed by employing the 
iteration method proposed in [3lj , which is based on the 
Lee- Suzuki similarity transformation [32| . 

The Viow-k given by the T-matrix equivalence ap- 
proach mentioned above is not hermitian. Therefore, an 
additional transformation is needed to make it hermi- 
tian. To this end, we resort to the hermitization proce- 
dure suggested in [3l|, which makes use of the Cholesky 
decomposition of symmetric positive definite matrices. 



A. Potential model Vi ow -k 

First, we outline the derivation of Viow-k [l3l. [13. flU- 
As pointed out in the Introduction, the repulsive core 
contained in Vnn is smoothed by integrating out the 
high-momentum modes of Vnn down to a cutoff momen- 
tum A. This integration is carried out with the require- 
ment that the deuteron binding energy and low-energy 
phase shifts of Vnn are preserved by Vj ow -k- This is 
achieved by the following T-matrix equivalence approach. 
We start from the half-on-shell T matrix for Vnn 

T(k',k,k 2 ) = V NN (k',k) + 

f°° 1 

+ V / q 2 dqV N N{k',q)- 2 ^T(q,k,k 2 ) , (I) 

J o k — q 

where V denotes the principal value and k, k! , and 
q stand for the relative momenta. The effective low- 
momentum T matrix is then defined by 



Tiow-k(p',P,p 2 ) = Mow-k(p',f>) + 



,A l 

-V / q 2 dqV\ ow - k (p' ,q)-z -T iovj _ k (q,p,p 2 

Jo p -q 



(2) 



where the intermediate state momentum q is integrated 
from to the momentum space cutoff A and (p',p) < A. 
The above T matrices are required to satisfy the condi- 
tion 

T(p',p,p 2 ) = T low _ k (p',p,p 2 ); (p',p)<A. (3) 

The above equations define the effective low- 
momentum interaction Vq ow -kj and it has been shown 
[lij that they are satisfied by the solution: 

VW-k - Q-Q' f Q+Q' I Q I Q-Q 1 f Q f Q [ Q+- 



(4) 

which is the well known Kuo-Lee-Ratcliff (KLR) folded- 
diagram expansion [l], originally designed for con- 
structing shell-model effective interactions. In Eq. (|U) 



B. Potential model N^LOW 

The general and fundamental reason why the Viow-k 
approach to nuclear structure physics works is that the 
dynamics which rules nuclear physics can be described in 
the framework of a low-energy EFT. This nuclear EFT 
is characterized by the symmetries of low-energy QCD, 
in particular, spontaneously broken chiral symmetry, and 
the degrees of freedom relevant for nuclear physics, nu- 
cleons and pions. The expansion based upon this EFT 
has become known as chiral perturbation theory (%PT), 
which is an expansion in terms of (Q/Mqcd) u where Q 
denotes the magnitude of a nucleon three-momentum or 
a pion four-momentum, and Mqcd is the QCD energy 
scale [J, .5]. For this expansion to converge at a proper 
rate, we have to have Q <§; ATqcd- To enforce this, chi- 
ral NN potentials are multiplied by a regulator function 
that suppresses the potential for nucleon momenta Q > A 
with A < A/qcd- Present chiral NN potentials [Til. Il2j 
typically apply values for A around 2.5 fin -1 . 

It is of course not an accident that the latter cutoff 
value is not too different from what is typically used for 
Viow-k; namely, A = 2.1 fm _1 . This fact stimulates an 
obvious question, namely: Is it possible to construct a 
chiral NN potential with A = 2.1 fru -1 ? Not surpris- 
ingly, the answer is in the affirmative. 

Thus, we have constructed a NN potential at N 3 LO of 
chiral perturbation theory that carries a sharp momen- 
tum cutoff at 2.1 fm _1 . We have dubbed this potential 
N 3 LOW. 

One advantage of this potential is that it is given in 
analytic form. The analytic expressions are the same as 
for the "hard" N 3 LO potential by Entem and Machleidt 
of 2003 [TO] and are given in Ref. [H]. Note that the 
procedure that needs to be followed to construct Viow-k 
from a free NN potential (cf. Sec. II. A) can be carried 
out only numerically. 

Another important issue are many-body forces. As 
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demonstrated in Ref. [301 ] , when applying a Vi ow -k m cer- 
tain nuclear many-body systems, the inclusion of three- 
body forces (3NF) may be crucial. For example, nuclear 
matter does not saturate without a 3NF when the two- 
nucleon force is represented by a low-momentum poten- 
tial. One great advantage of xPT is that it generates 
nuclear two- and many-body forces on an equal footing 
(for an overview of this aspect, see Ref. [3J|). Most 
interaction vertices that appear in the 3NF and in the 
four-nucleon force (4NF) also occur in the two-nucleon 
force (2NF). The parameters carried by these vertices are 
fixed in the construction of the chiral 2NF. Consistency 
requires that for the same vertices the same parameter 
values are used in the 2NF, 3NF, 4NF, .... If the 2NF 
is analytic, these parameters are known, and there is no 
problem with their consistent proliferation to the many- 
body force terms. However, if a potential exists only 
in numeric form, then those parameters are not explic- 
itly known and the parameters to be used in the 3NF, 
4NF, ...must be based upon 'educated guesses'. The 
firm consistency between two- and many-body forces is 
lost. Moreover, if a potential is given only in numeric 
form, one does not known what order of xPT it belongs 
to. Thus, it is also not clear which orders of 3NF and 
4NF to include to be consistent with the order of the 
2NF. 
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FIG. 2: Same as Fig. \T\ 
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FIG. 1: Phase parameters of neutron-proton scattering up to 
200 MeV laboratory energy for partial waves with total an- 
gular momentum J < 2. The solid lines show the predictions 
by the new low-momentum N 3 LOW potential. The solid dots 
and open circles represent the Nijmegen multienergy np phase 
shift analysis [5l[ and the GWU/VPI single-energy np anal- 
ysis SM99 [52]|. respectively. 



Our newly constructed N 3 LOW reproduces accurately 
the empirical deuteron binding energy, the experimen- 
tal low-energy scattering parameters, and the empirical 
phase-shifts of NN scattering up to at least 200 MeV 
laboratory energy, see Figs. 11121 More details about this 
potential will be published elsewhere 35]. It is the main 
purpose of this paper to test if the perturbative properties 
of N 3 LOW when applied in microscopic nuclear structure 
are as good as the ones of typical Vjow-k potentials. 



III. DERIVATION OF THE SHELL-MODEL 
EFFECTIVE HAMILTONIAN 

In the framework of the shell model, an auxiliary one- 
body potential U is introduced in order to break up the 
nuclear hamiltonian as the sum of a one-body compo- 
nent Hq, which describes the independent motion of the 
nucleons, and a residual interaction Hi: 



^ 2m ^ 



(T + U) + (V - U) = H + H x 



(5) 



Once Hq has been introduced, a reduced model space 
is defined in terms of a finite subset of H 's eigenvectors. 
The diagonalization of the many-body hamiltonian (O 
in an infinite Hilbert space, that is obviously unfeasible, 
is then reduced to the solution of an eigenvalue problem 
for an effective hamiltonian H Q g in a finite space. 
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In this paper we derive H e g by way of the time- 
dependent perturbation theory [J, 0]. Namely, H e g is 
expressed through the KLR folded-diagram expansion in 
terms of the vertex function Q-box, which is composed 
of irreducible valence-linked diagrams. We take the Q- 
box to be composed of one- and two-body Goldstone di- 
agrams through third order in V [36j| . Once the Q-box 
has been calculated at this perturbative order, the series 
of the folded diagrams is summed up to all orders using 
the Lee-Suzuki iteration method [32[. 

The hamiltonian H c q contains one-body contributions, 
which represent the effective SP energies. In realistic 
shell-model calculations it is customary to use a subtrac- 
tion procedure so that only the two-body terms of 
H e g are retained - the effective interaction V c s - and the 
SP energies are taken from the experimental data. 

In this work, we have followed a different approach 
and employed the theoretical SP energies obtained from 
the calculation of H c s . This allows to make a consistent 
study of the perturbative properties of the input poten- 
tial V (Vi ow _ k or N 3 LOW) in the shell-model approach. 

In this regard, it is worth to point out that, owing to 
the presence of the — U term in Hi, [/-insertion diagrams 
arise in the Q-box. In our calculation we use as auxiliary 
potential the harmonic oscillator (HO), U = ^mcj 2 r 2 +A, 
and take into account only the first-order [/-insertion 
that appears in the collection of the one-body contribu- 
tions, which is the dominant one 38]. 



IV. RESULTS 

We have performed shell-model calculations for even- 
mass oxygen isotopes beyond the doubly-closed core 16 O. 
Oxygen isotopes constitute a quite interesting nuclear 
chain both theoretically and experimentally, and they 
have been considered a testing ground for realistic shell- 
model calculations since the pioneering work by Kuo and 
Brown [3f|. Our calculations have been carried out by 
using the Oslo shell-model code (4^| 

Two -ffeff's have been obtained using, respectively, a 
Mow-k with a cutoff momentum A = 2.1 fm _1 derived 
numerically from the 'hard' N 3 LO chiral NN potential 
(llj . and the new N 3 LOW potential described in Sec. 
II. B. The Coulomb force between proton-proton interme- 
diate states has been explicitly taken into account. For 
the HO parameter fuv we have used the value of 14 MeV, 
as obtained from the expression hu> — 45A -1 / 3 — 25 A~ 2 ^ 3 
[4l| for A — 16. The value of A has been chosen so 
that the self-energ y an d the [/-insertion diagrams almost 
cancel each other [38j]. The numerical value is -54 MeV 
yielding an unperturbed energy of the Osld shell equal to 
-5 MeV, which is not far the experimental value of the 
ground state (g.s.) energy of 17 relative to 16 (-4.144 
MeV [HI). However, it is worth to point out that, when 
including only first-order [/-insertions, our H e ff matrix 
elements do not depend on the choice of A. 

The perturbative behavior of the input potential V 



rules the convergence rate of the diagrammatic series for 
H e ff which has to deal with the convergence of both the 
order-by-order perturbative expansion and the sum over 
the intermediate states in the Goldstone diagrams. In 
this context, we have found it interesting to study and 
compare the convergence properties of H c g derived from 
both the Viow-k and the N 3 LOW potentials. 

TABLE I: Ground state energies (in MeV) of 18 relative to 
16 O calculated with H e g derived from the VW-k of the "hard" 
N 3 LO potential as a function of the maximum number N max 
of the HO quanta (see text for details). Results obtained at 
second and third order in perturbation theory are reported. 

TV^ 4 (3 8 10 12 14 W~ 

2nd -8.191 -10.615 -12.748 -14.318 -15.037 -15.142 -15.168 
3rd -6.617 -8.987 -11.344 -13.277 -14.294 -14.487 -14.523 



TABLE II: Same as Table |U but with H cS derived from 
N 3 LOW. 



iV max 4 6 8 10 12 14 16 
2nd -4.806 -6.789 -9.085 -11.759 -13.671 -14.108 -14.162 
3rd -3.547 -5.366 -7.683 -10.789 -13.339 -13.992 -14.049 



In Table U we present the Vi ow _k g.s. energies of 18 
relative to 16 O obtained at second- and third-order per- 
turbative expansion of H e g . In both cases the en- 
ergies are reported as a function of the maximum al- 
lowed excitation energy of the intermediate states, ex- 
pressed in terms of the oscillator quanta N max . It can 
be seen that the g.s. energy is practically convergent 
at N max — 12. As regards the order-by-order conver- 
gence, it may be considered quite satisfactory, the dif- 
ference between second- and third-order results being 
around 4 % for N max = 16, which is the largest /V max 
we have employed. Similar results are obtained using the 
N 3 LOW potential, as shown in Table |TTJ In this case 
for N max = 16 the difference between second- and third- 
order results is less than 1 %. 

TABLE III: Same as Table HI but with experimental SP ener- 
gies. 

4 6 8 10 12 14 16" 

2nd -12.201 -12.000 -11.847 -11.771 -11.750 -11.748 -11.749 
3rd -12.373 -12.328 -12.286 -12.281 -12.290 -12.294 -12.296 



TABLE IV: Same as Table HH but with experimental SP en- 
ergies. 



-^Vmax 


4 


6 


8 


10 


12 


14 


16 


2nd 


-12.386 


-12.202 


-12.072 


-12.002 


-11.973 


-11.968 


-11.967 


3rd 


-12.663 


-12.637 


-12.612 


-12.656 


-12.722 


-12.737 


-12.739 



It is now worth to comment on the dimension of the 
intermediate-state space. We have found that, when in- 
creasing -/V max , the matrix elements of the two-body ef- 
fective interaction V e g (TBME) are much more stable 
than the calculated SP energies. In this connection, an 
inspection of Tables IIIII and IIVI where we present re- 
sults obtained using experimental SP energies, evidences 
the rapid convergence of the g.s. energy with N max . 
This supports the choice of moderately large interme- 
diate state spaces in standard realistic shell-model cal- 
culations, where in general experimental SP energies are 
employed. 

A comparison between the results obtained from the 
Mow-k and the N 3 LOW potentials (see Tables Q] and [Til 
respectively) shows that at third order, with a sufficiently 
large number of intermediate states, the two interactions 
predict very close values for the relative binding energy 
of 18 0. Note that the experimental value is -12.187 MeV 
|42| . From now on we refer to calculations with N max — 
16 and Q-box diagrams up to third order in perturbation 
theory. 



TABLE V: r = 0+ TBME of H cS (in MeV) obtained from 
the Viow-k of the "hard" N 3 LO potential and N 3 LOW at third 
order in perturbation theory with iV max = 16. They are com- 
pared with the USD A TBME of Brown and Richter [H|. 



configuration 


Fiow-k 


N 3 LOW 


USDA 


(0d5/2) 2 (0d5/2) 2 
(0d5/2) 2 (0d3/2) 2 
(0d5/2) 2 (lsl/2) 2 
(0d3/2) 2 (0d3/2) 2 
(0d3/2) 2 (lsl/2) 2 
(lsl/2) 2 (lsl/2) 2 


-2.435 


-2.689 


-2.480 


-3.464 


-3.741 


-3.569 


-1.269 


-1.337 


-1.157 


-0.845 


-1.147 


-1.505 


-0.862 


-0.967 


-0.983 


-2.385 


-2.563 


-1.846 



USDA hamiltonian, see Ref. [43| for details. It is inter- 
esting to note that the TBME of USDA are closer to ours 
than those of the original USD [HI, [4{| , which was based 
on a smaller set of data. Regarding the SP energies, we 
see that the USDA absolute energy of the d 5 / 2 state is in 
good agreement with the experimental g.s. state energy 
of 17 0, while both our calculations give about one MeV 
more binding. On the other hand, our calculated and the 
USDA excitation energies of the Si/ 2 state come close to 
each other and do not differ significantly from the exper- 
imental value. As regards the relative energy of the d 3 / 2 
state, the USDA value is larger than the experimental 
one, in line with our findings. 
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FIG. 3: (Color online) Experimental and calculated ground- 
state energy per valence neutron for even oxygen isotopes 
from A — 18 to 28. iV va i is the number of valence neutrons. 
The calculated values are obtained with H B g derived from 
the Fiow-k of the "hard" N 3 LO potential at third order in 
perturbation theory, with 7V max = 16. 



TABLE VI: Calculated SP relative energies of H ef! (in MeV) 
obtained from the Fi ow -k of the "hard" N 3 LO potential and 
N 3 LOW at third order in perturbation theory with iV max = 
16. They are compared with the USDA SP energies of Brown 
and Richter [431 ] and the experimental ones. The values in 
parenthesis are the absolute SP energies. 



orbital 


Fiow-k 


N 3 LOW 


USDA 


Expt 


vQdh/2 


0.0 (-5.425) 


0.0 (-4.909) 


0.0 (-3.944) 


0.0 (-4.144) 


i^0d3/2 


7.323 


7.117 


5.924 


5.085 


vlsl/2 


1.257 


0.818 


0.883 


0.871 



In Table [V] and EO we present the J = 0+ TBME 
and the SP energies obtained from the Fj ow -k and the 
N 3 LOW potentials. It turns out that the results from 
the two potentials are quite similar, the largest difference 
regarding the absolute energy of the d 5 / 2 SP level. In the 
last column of Table fVl and I VII we also show the = + 
TBME and SP energies obtained by Brown and Richter 
[43j from a least-squares fit of a large set of energy data 
for the sd-shell nuclei. The reported values refer to the 
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FIG. 4: (Color online) Same as Fig. [3] but with H c g derived 
from N 3 LOW. 
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FIG. 5: (Color online) Same as Fig. [3] but with SP energies 
shifted upward by 1.1 MeV. 
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mental 18 O binding energy, then we obtain a very good 
agreement, as shown in Figs. [5] and [51 

This result is very interesting: our calculations fail to 
predict that 26 O and 28 O are unbound to 2-neutron decay 
as indicated by experimental studies (46l. |47|. However, 
when the SP energies are shifted upward, then we obtain 
the correct binding properties, as it may be seen from 
Figs. and [5J In some recent papers [4g, |49( it has 
been argued that with a realistic shell-model interaction 
it is not possible to reproduce the correct behavior of the 
binding energy of oxygen isotopes. In Refs. [4!j[5(| this 
has been ascribed to the lack of both genuine and effective 
three-body correlations. Our results suggest that a major 
factor is the lack of a real three-body force. In fact, as 
discussed above, we have obtained a good description 
of the binding energy of oxygen isotopes by modifying 
the SP spectrum so as to reproduce the g.s. energy of 
two- valence nucleon 18 O. The last quantity is obviously 
not affected by effective three-body correlations, but it is 
sensitive to genuine three-body forces. 

In Figs. [7] and [5] we compare the experimental and 
calculated excitation energies of the first 2 + states for 
the I4ow-k and the N 3 LOW potential, respectively. In 
both cases the experimental behavior as a function of A 
is well reproduced. 



7.0 
6.0 
5.0 
4.0 
3.0 
2.0 
1.0 
0.0 







♦ Expt. 




• Theor. 





16 



20 22 
A 



24 



26 



FIG. 6: (Color online) Same as Fig. [4] but with SP energies 
shifted upward by 1 MeV. 



It is worth recalling that -ff c ff , which has been derived 
in the framework of the KLR approach, is defined for the 
two valence- nucleon problem. It is of interest, however, 
to test the theoretical SP energies, as given in Table PVTl 
and the TBME in systems with many valence nucleons. 
To this end, we have calculated for even oxygen isotopes 
the g.s. energies relative to 16 core up to A=28, and 
the excitation energies of the first 2 + state up to A=24. 
In Figs. [3] and 0] the experimental g.s. energies per va- 
lence nucleon are compared with the calculated ones for 
the V[ ow -k and the N 3 LOW potential, respectively. The 
predictions from both potentials yield binding energies 
which are larger than the experimental ones, with the 
difference remaining almost constant when increasing the 
number of valence neutrons N va \. In fact, if we shift the 
VIow-k and N 3 LOW SP spectra upward by 1.1 MeV and 
1 MeV, respectively, such as to reproduce the experi- 



FIG. 7: (Color online) Experimental and calculated excitation 
energy of the first 2 + state for oxygen isotopes. The calculated 
values are obtained with H c g derived from the Viow-k of the 
"hard" N 3 LO potential. Experimental data are taken from 
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FIG. 8: (Color online) Same as in Fig. [7] but with H e a 
derived from N 3 LOW. 
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V. SUMMARY AND CONCLUSIONS 

In this paper, we have performed shell-model calcula- 
tions employing an effective hamiltonian obtained from a 
new low-momentum potential, dubbed N 3 LOW, which is 
derived in the framework of chiral perturbation theory at 
next-to-next-to-next-to-leading order with a sharp cutoff 
at 2.1 far 1 . 

We have studied the convergence properties of this po- 
tential by calculating the ground-state energy of 18 O, 
and have compared the results with those obtained using 
the Viow-k derived from the "hard" N 3 LO potential of 
Entem and Machleidt [ll|. It has turned out that the 
two low-momentum potentials show the same perturba- 
tive behavior. We have also performed calculations for 
the heavier oxygen isotopes, obtaining binding energies 
which are very similar for the two potentials, but larger 
than the experimental ones. However, we have found out 
that by introducing an empirical shift of the SP spectrum 
the ground-state energies are well reproduced up to 24 0, 
while 26 O and 28 O are unbound to two-neutron decay, in 
agreement with the experimental findings. 



The main purpose of this work was to test the new 
chiral low-momentum NN potential N 3 LOW. Unlike 
Viow-k, this potential has the desirable feature that it 
is given in analytic form. We have shown here that it 
is suitable to be applied perturbatively in microscopic 
nuclear structure calculations and yields results which 
come quite close to those obtained from the Ui ow -k de- 
rived from the "hard" N 3 LO potential. 

After this first successful application, there is now mo- 
tivation to pursue further nuclear structure calculations 
with this new low-momentum potential. 
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